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IV 


1.  SUMMARY 


Yield  estimation  of  small  explosions  at  local  distances  represents  a  challenge  for  the 
nuclear  explosion  monitoring  community.  We  have  examined  the  feasibility  of  using  short- 
period  surface  wave  magnitudes,  called  MRg,  to  estimate  explosion  yields  at  local  distances  (< 
100  km).  We  have  modified  the  Russell  (2006)  Ms  formula,  which  was  derived  for  periods  of  8- 
25  s  for  distances  beyond  50  km,  for  application  at  local  distances  and  1  s  period.  We  have 
studied  short-period  surface  wave  attenuation  in  diverse  lithologies  in  order  to  incorporate  an 
attenuation  term  in  the  magnitude  scale,  which  is  suitable  for  Rg.  We  have  also  incorporated 
multiple  excitation  corrections  for  Rg  based  on  the  near  source  seismic  velocities,  which  greatly 
affect  the  source  region  amplitudes  for  Rg.  It  is  important  to  note  that  in  the  formula  the 
excitation  is  estimated  from  the  measured  Rg  group  velocity.  We  have  also  derived  a  new 
Butterworth  filter  cutout  definition  for  filtering  Rg  near  1  s  period  at  distances  between  2-100 
km.  We  used  the  new  formula  to  estimate  MRg  for  39  small  (37  <  Y  <  12,270  kg  TNT 
equivalent)  and  shallow  (<120  m)  explosions  detonated  in  North  America  in  lithologies  ranging 
from  alluvium  to  granite.  Regressions  of  the  magnitudes  with  yield  result  in  the  equation  MRg 
=logio(Y)-3.03  for  chemical  explosions  with  Y  <  12,270  kg.  An  F  factor  with  95%  confidence 
was  determined  to  be  2.25,  giving  lower  and  upper  bounds  on  the  yield  estimates  of  Y/2.25  and 
Y*2.25,  respectively.  We  applied  the  relationship  MRg  =logio(Y)-3.33  (assuming  factor  of  2 
equivalence  between  chemical  and  nuclear)  to  nuclear  explosions  detonated  at  the  Degelen  and 
Shagan,  Kazakhstan,  test  sites.  The  estimated  yields  based  on  Rg  magnitudes  were  often  within 
20%  of  the  true  yield  and  had  smaller  F  factor  than  the  estimated  yields  for  US  chemical 
explosions. 
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2.  INTRODUCTION 


The  accurate  estimation  of  yields  from  small  chemical  and  nuclear  explosions  represents 
a  difficult  and  continuing  problem  for  the  Nuclear  Explosion  Monitoring  (NEM)  community.  P- 
wave  spectra  (Murphy,  1996),  Pn  amplitudes  (Vergino  and  Mensing,  1990),  Lg  amplitudes 
(Nuttli,  1986),  intermediate-period  surface  wave  amplitudes  (Stevens  and  Murphy,  2001),  or 
coda  techniques  (Mayeda  et  al.,  2003)  are  some  of  the  methods  that  have  been  proposed  for 
estimating  the  yields  of  small  explosions  recorded  at  regional  distances.  However,  the 
uncertainty  on  seismic  yield  estimates  can  be  large,  and  reducing  the  uncertainty  may  require  a 
priori  information  about  source  media,  knowledge  of  the  emplacement  depth,  and  calibrations 
for  path  and  site  effects.  There  is  no  general  consensus  as  to  which  method  for  seismic  yield 
estimation  works  best  for  all  nuclear  test  sites. 

For  monitoring  at  local  and  near-regional  distances,  we  propose  a  methodology  for  yield 
estimation  based  on  magnitudes  of  short-period,  fundamental  mode  Rayleigh  waves  ( Rg ).  At 
local  distances,  Rg  can  be  the  largest  amplitude  seismic  arrival  observed  from  shallow 
explosions,  mining  explosions,  and  shallow  earthquakes.  Adushkin  (2001)  demonstrated  the 
ability  to  use  Rg  amplitudes  at  near-regional  distances  to  provide  accurate  seismically-estimated 
yields.  He  corrected  Rg  amplitudes  from  Semipalatinsk  underground  explosions  for  attenuation, 
geometric  spreading,  and  station-specific  effects.  Using  the  corrected  Rg  amplitudes,  he 
estimated  explosive  yield  within  20%  of  the  true  yields  for  the  explosions. 

In  this  paper,  we  develop  a  new  formula  for  short-period  (~1  s)  surface  wave  magnitudes, 
called  Ms#,  using  the  methods  of  Russell  (2006).  Attenuation  coefficients  for  short-period 
surface  waves  have  been  derived  from  diverse  tectonic  settings  to  calibrate  the  new  M«A,  formula. 
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We  demonstrate  the  utility  of  the  formula  at  estimating  magnitudes  for  small  explosions,  and 
then  correlate  the  magnitudes  to  yield  with  uncertainty  estimates. 


3.  MAGNITUDE  FORMULA  DEVELOPMENT 


Introduction 

While  many  surface  wave  magnitude  scales  have  been  developed  during  the  past  century, 
none  were  developed  specifically  for  Rg.  Russell  (2006)  developed  a  magnitude  formula  that 
could  effectively  measure  surface-wave  magnitudes  at  local,  regional  and  teleseismic  distances, 
at  variable  periods  between  8  and  25  seconds.  With  suitable  modifications,  we  will  extend  this 
method  to  be  applicable  to  1  sec  filtered  Rg  phases.  The  (2006)  magnitude  equation  is: 


where  ah  is  the  amplitude  of  the  Butterworth-filtered  surface  waves  (zero-to-peak  in 


nanometers),  A  is  the  distance  in  degrees,  and  fc  is  the  filter  cutoff  frequency  of  3rd  order,  zero 


phase  Butterworth  filters  with  corner  frequencies  of  l/T-fc,  l/T+fc,  respectively. 

For  periods  8<T<25,  the  equation  is  corrected  to  T=  20  s,  accounting  for  frequency- 
dependent  source  effects,  attenuation,  and  dispersion.  The  constant,  C,  is  determined  to  be  0.43, 
which  scales  the  equation  at  T=20  s  to  be  equivalent  to  Rezapour  and  Pearce  (1998). 

In  the  next  few  paragraphs,  we  discuss  how  Eq.  1  can  be  modified  for  short-period 
Rayleigh  waves  magnitudes,  which  we  refer  to  as  MRg.  This  will  include  a  change  in  the  filter 
order  which  results  in  a  new  constant  in  Eq.  1,  new  filter  constraints  to  determine  fc,  a  new 
attenuation  correction  based  on  observed  Rg  data  in  diverse  geologies,  and  choices  for  the 
excitation  correction  based  on  near-source  medium  properties.  These  changes  will  also  scale 
corrections  to  standard  20  sec  Ms  as  in  Eq.  1. 
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Filter  Order.  The  first  modification  to  the  Russell  (2006)  fonnula  for  MRg  is  to  use  a  2nd  order 
Butterworth  filter  rather  than  3rd  order.  This  decision  was  based  on  multiple  comparisons  of 
filtered  Rg  waveforms  using  2nd  and  3rd  order  filters,  and  preferring  2nd  order  filtering  in  the 
identification  and  processing  of  the  waveforms  at  local  distances.  The  constant  (C=0.43)  in  Eq.  1 
is  changed  to  C=0.46  as  a  result  of  applying  a  2nd  order  filter,  due  to  the  constant  being 
dependent  on  the  Butterworth  gain  term  f?n,  which  is  a  function  of  the  Butterworth  filter  order 
(Russell,  2006). 

Filter  Definition.  To  estimate  the  time-domain  amplitude  of  Rg,  we  filter  the  time  series  with  a 
2nd  order,  zero-phase  Butterworth  band  pass  filter  with  corner  frequencies  1/T-fc,  1/T+fc,  where 


0.1 

7VA 


(2) 


The  constant,  0.1,  in  Eq.  2  is  based  on  the  dispersive  properties  of  short-period  Rayleigh  waves, 
which  often  exhibit  strong  dispersion  in  the  T=1  s  range.  For  8-25  s  Rayleigh  waves  (e.g., 
Russell,  2006),  the  constant  was  0.6.  We  note  that  Zeiler  and  Velasco  (2009)  used  a  constant  of 
0.2  in  Eq.  2  to  extend  Eq.  1  to  short-period  Rg  for  an  explosion  dataset  in  Arizona. 

Rg  Attenuation.  We  estimated  the  Rg  attenuation  coefficient,  y,  for  explosion  datasets  in 
southern  Arizona  (igneous  rocks),  Colorado  Plateau  (sedimentary),  central  Alaska 
(metamorphic),  New  England  (granite),  Albuquerque,  New  Mexico  (alluvium  and  sedimentary), 
Yucca  Flat  in  Nevada  (alluvium),  and  Wyoming  (sedimentary).  The  period-dependent 
attenuation  coefficients  for  all  datasets  (Figure  1)  show  that  for  most  of  the  explosion  datasets, 
intrinsic  attenuation  of  Rg  results  in  an  average  y=0.05±0.03  /km  at  T=1  s.  We  note  that  this 
value  will  not  be  sufficient  in  an  area  of  strong  Rg  scattering,  such  as  near  faults  (see  Figure  1 
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solid  data  points)  or  rapid  topographic  changes.  This  is  not  a  significant  problem  given  that 
identifiable  Rg  will  likely  not  propagate  to  many  stations  beyond  such  fault  zones. 

Figure  1  also  shows  period-dependent  attenuation  correction  Eq.  1  of  the  form: 

Ba,t(To/T)P  (3a) 

where  Battand  p  were  0.0031  and  1.8  for  8-25  s  surface  waves.  We  find  that  to  better  match  the 
Rg  attenuation  coefficients  at  T=0.1  to  2  s,  the  exponent  (p)  value  should  be  increased  to  2.19 
which  corresponds  to  a  y  =  0.05  /km  at  T=1  s.  We  initially  used  this  value  in  the  development  of 
MRg,  but  noted  that  for  our  small  chemical  explosion  dataset,  if  Rg  was  observed  beyond  ~25 
km,  it  was  usually  in  a  region  with  lower  attenuation  coefficients.  As  a  result,  the  MRg  at 
distances  greater  than  20  km  would  be  overcorrected  for  attenuation  using  a  y  =  0.05  /km.  Thus 
we  chose  to  weight  the  attenuation  correction  in  MRg  to  a  lower  y  =  0.02  /km  (or  an  equivalent 
exponent  of  p=1.96).  Returning  to  Eq.  3a,  the  attenuation  tenn  for  our  proposed  MRg  formula  at 
T=  1  sec  is: 

0.003 1(20/T)L96A  ,  (3b) 
or  1.1  A.  (3c) 

Rg  Excitation  Correction.  Eq.  1  includes  a  source  excitation  correction  of  — ii.66  Log  in  the 

estimation  of  surface  wave  magnitudes  at  T=8-25  s.  The  correction  is  used  to  reduce  short-period 
Rayleigh-wave  amplitudes,  which  are  enhanced  for  shallow  sources — particularly  explosions — 
in  order  to  approximate  a  magnitude  that  would  be  measured  for  the  same  event  near  the 
reference  period  of  T=20  s.  Figure  2  shows  the  theoretical  magnitude  correction  required  to 
extend  the  excitation  correction  to  T=1  s  Rg  for  five  different  crustal  models.  The  corrections  are 
based  on  explosion  synthetics  at  1  km  distance  for  a  source  at  0.01  km  depth.  These  corrections 
are  normalized  at  T=20  s  for  scaling  purposes.  Also  shown  is  the  excitation  correction  from  Eq. 
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1,  which  matches  the  theoretical  magnitude  correction  for  most  models  in  the  8-25  s  range.  We 
note  that  while  these  corrections  are  for  a  synthetic  explosion  at  0.01  km  depth,  we  find  similar 
corrections  up  to  about  200  m  depth,  except  for  the  alluvium  correction,  which  is  slightly 
reduced. 

Figure  2  suggests  that  a  single  excitation  correction  may  not  provide  optimal  results  for 
estimating  MRg  magnitudes  in  differing  geological  environments.  The  magnitude  corrections 
range  from  -0.6  m.u  for  the  global  average  AK135  model  (Kennett  et  ah,  1995)  to  -1.4  m.u.  for  a 
western  United  States  (WUS)  crustal  model  (Herrmann  et  al.,  2011).  If  the  upper  km  of  the 
WUS  model  is  replaced  with  velocities  representing  dry  alluvium,  the  resulting  excitation 
correction  increases  to  approximately  -2.6.  The  nulls  in  the  spectrum  for  the  alluvium  model 
show  the  sensitivity  of  the  synthetics  to  the  low  velocities.  The  Eq.  1  correction  term  is  -0.9  at  T 
=  1  s  which  is  similar  to  the  theoretical  correction  based  on  velocity  models  for  the  Korean 
Peninsula  (Cho  et  al.,  2007)  and  the  central  United  States  (Hernnann  et  al.,  2011). 

The  velocities  of  the  layer  where  the  explosion  is  synthesized  primarily  controls  the 
magnitude  correction.  For  example,  the  Vp  and  Vs  for  the  upper  layer  of  the  western  United  States 
model  are  3.4  km/s  and  2  km/s,  respectively.  Conversely,  for  the  AK135  model,  the  VP  and  Vs  are 
5.8  km/s  and  3.46  km/s,  respectively.  For  the  WUS+dry  alluvium  model,  the  Vp  and  Vs  for  the 
upper  layer  are  2.0  and  0.8  km/s,  respectively.  Based  on  these  velocities,  we  propose  the 
following  excitation  correction,  E(T=1  s),  for  the  MRg  fonnula: 

-2.5  UK(f  <  0.5  km/s 
E(T=1  s)  =  -1.5  0.0  <  ^  2  km/s  ,  (4) 

—  1  UK(S  >  2  km/s 

where  URg  is  the  measured  group  velocity  of  T=1  s  Rg  near  the  test  site.  This  highlights  the 
importance  of  emplacement  properties  on  the  generation  of  Rg  and  provides  for  a  larger 
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magnitude  correction,  relative  to  T=20  s,  for  alluvium  and  “soft”  rock  test  sites  than  “hard”  rock 
test  sites.  It  requires  an  approximate  origin  time  and  distance  in  order  to  process  for  MRg.  It 
should  be  emphasized  that  the  corrections  in  Eqn  (4)  are  based  on  the  assumption  that  at  local 
source/receiver  distances,  the  source  geology  is  the  same  as  for  the  propagation  path.  For  more 
complex  media  where  the  source  group  velocities  differ  from  the  propagation  velocities, 
knowledge  of  the  source  group  velocity  from  independent  information  will  be  necessary  to 
correct  Eqn.  (4). 

Final  Mr§  Formula 

Based  on  the  described  changes  to  the  Russell  (2006)  fonnula,  we  introduce: 

Mro  =  log(aRg)  +  Vi  log  (sin  (A))  -  log {fc]  +  1.1A  +  E(T)  -  0.46  (5) 

where  A  is  distance  in  degrees  between  0.02°  and  0.90°  (2  and  100  km).  The  current  distance 
limitations  are  constrained  by  the  data  that  are  discussed  in  this  paper.  Waveforms  should  be 
corrected  for  instrument  response  and  converted  to  displacement  (in  nanometers)  and  filtered 
between  1  -fc,  1  +fc  using  a  zero-phase,  2nd  order  Butterworth  bandpass  filter  where  /■  is 
determined  by  Eq.  2.  Zero-to-Peak  amplitudes  of  Rg  (aRg)  in  nanometers  are  measured,  and  the 
time  of  the  peak  amplitude  recorded  and  converted  to  group  velocity  (URg)  in  order  to  pick  an 
excitation  correction,  E(T),  defined  in  Eq.  4.  The  processing  can  be  automated  using  pre-selected 
group  velocity  windows  for  Rg  or  manually  by  an  analyst  picking  the  peak  amplitudes  for  Rg. 


4.  MRg  APPLICATION 

Explosions.  We  have  applied  the  MRg  fonnula  to  several  explosion  datasets  recorded  at  local 
distances.  Figure  3a  shows  the  Butterworth-filtered  waveforms  for  a  101.4  kg  (TNT  equivalent) 


Approved  for  public  release;  distribution  is  unlimited. 

7 


explosion  detonated  in  Barre  granite  and  recorded  from  2  to  24  km.  The  Butterworth  filtered 
waveforms  (Eq.  2)  highlight  the  T=1  s  Rg.  Group  velocities  of  the  peak  amplitude  of  the  filtered 
signals  are  between  2.2  and  2.8  km/s  and  increase  with  distance  (Figure  3b).  The  mean  Urs  is  2.6 
km/s,  thus  the  processing  automatically  chose  an  excitation  correction  of  E(T)=-1  for  Ur<,  >  2.0 
km/s,  which  is  typical  of  “hard”  rock  emplacement  media.  The  MRg  estimates  using  Eq.  5  are 
shown  in  Figure  3c.  The  network  average  Mr^,  which  is  estimated  after  discarding  ±5%  of  high 
and  low  estimates  in  the  population,  is  -1.05  with  of  ±  one  standard  deviation  of  0.07  m.u.  To 
ensure  that  we  are  not  biasing  these  results  by  using  instrument-corrected  4.5  Hz  geophone  data, 
we  calculated  the  magnitude  using  15  broadband  seismometers  and  obtained  MRg  =  -1.06±0.10. 

We  also  estimated  MRg  for  39  explosions  recorded  at  distances  between  2  and  50  km. 
The  explosions  were  part  of  multiple  active  source  phenomenology  experiments,  including  the 
Arizona  Source  Phenomenology  Experiment  (AZ-SPE),  Bighorn  Active  Seismic  Experiment 
(BASE),  the  HUMMING  ALBATROS  Experiment  (HAL),  the  Fracture  Decouple  Experiment 
(FDE),  the  Frozen  Rock  Experiment  (FRE),  HUMBLE  REDWOOD  (HR)  I  and  II,  and  the  New 
England  Damage  Experiment  (NEDE).  The  explosion  yields  for  the  40  events  range  from  37  kg 
to  12,270  kg  TNT  equivalent.  All  of  the  explosions  were  detonated  below  ground;  however,  the 
explosions  include  both  fully-  and  partially-confined  shots.  All  of  the  explosions  were  conducted 
at  shallow  depths  of  burial  (DOBs)  less  than  120  m;  however,  the  approximate  scaled  DOBs 
ranged  from  43  to  271  m/kt1/3. 

Yield-MRg  Relationship.  Figures  4  and  5  compare  the  estimated  MRg  to  TNT  equivalent  yield. 
In  Figure  4a, b  the  results  for  each  experiment  are  plotted  as  a  different  color  or  symbol.  To  show 
the  effect  of  the  excitation  correction,  we  include  results  (Figure  4a)  for  MRg  with  three  possible 
excitation  corrections  as  defined  in  Eq.  4,  and  we  present  the  MRg  results  (Figure  4b)  using  a 
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single  correction  (-1)  for  all  events.  The  use  of  multiple  excitation  corrections  significantly 
decreases  the  variance  of  the  magnitude-yields  than  the  use  of  a  single  correction.  The  observed 
magnitudes  in  Figure  4b  confirm  the  theoretical  results  shown  in  Figure  2,  which  show  that  the 
slower  velocity  rocks  (alluvium,  sedimentary)  excite  larger  amplitude  T=1  s  Rg  than  a  similar 
yield  explosion  in  faster  media  (e.g.,  granite).  We  plan  to  continue  to  develop  this  formula  and 
the  excitation  corrections  as  more  data  become  available,  especially  shots  in  alluvium.  We 
regressed  log  yield  (Y)  versus  MRg  in  Figure  5  to  determine  the  relationship: 

MRg  =  logio(Y)-3.03  (Chemical)  (6) 

for  37  <  Y  <12270  kg  (TNT  equivalent).  The  slope  from  the  analysis  was  0.99,  thus  statistically 
representing  a  slope  of  unity  between  MRg  and  log  Y.  Slopes  for  historical  Ms-log  yield 
relationships  typically  fall  somewhere  between  0.75  to  1.2  (Murphy,  1977;  Bache,  1982; 
Evemden  and  Marsh,  1987).  Slopes  of  1  have  been  reported  by  Stevens  and  Murphy  (2001)  for 
the  Shagan  Test  Site  in  Kazakhstan  and  Bache  (1982)  for  multiple  datasets.  The  95%  confidence 
intervals  were  also  determined  for  the  data  and  plotted  in  Figure  5.  Note  that  the  MRg  /  yield 
formula  is  expressed  in  kilograms  of  yield;  to  change  to  kilotons  add  6  to  the  right  hand  side  of 
Eqn  (6). 

Another  way  of  determining  approximate  confidence  bounds  on  yields  of  explosions  is  to 
estimate  F  factors.  For  95%  confidence  interval,  the  F  factor  is  defined  as; 

F=10(1.96a/a)5  (y) 

where  a  is  the  standard  deviation  of  the  magnitudes  about  the  regression  line  with  slope  (a)  =  1 . 
The  regression  in  Figure  5  has  a  =  0.18  resulting  in  an  F=  2.25.  The  confidence  bounds  on  the 
yield  are  defined  as: 

Ylower  =Y/F  (8a) 
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Yupper=YxF  (8b) 

The  yield-MRg  relationship  (Eq.  6)  was  developed  for  a  series  of  chemical  explosions. 
The  source  function  for  the  chemical  explosion  is  equivalent  to  that  of  a  nuclear  explosion  with 
twice  the  yield  (Denny  et  al„  1995).  Thus  for  nuclear  explosions,  Eq.  6  would  be  modified  by 
subtracting  0.3  from  the  constant: 

MRg=logio(Y)-3.33  (Nuclear)  (9) 

Thus  for  a  1  kiloton  nuclear  explosion,  we  would  expect  an  MRg=2.7.  This  is  at  least  0.5  m.u. 
higher  than  most  historical  yield-  Ms  relationships  would  predict  (e.g.,  Stevens  and  Murphy, 

2001). 

Validation  of  the  Yield-MRg  Formula 

Kazakhstan.  We  have  validated  the  new  yield-MRg  fonnula  on  additional  datasets  that  were  not 
included  in  the  development  processes.  Validation  included  estimating  magnitudes  for  8 
Degelen,  Kazakhstan,  Test  Site  (DTS)  and  10  Shagan,  Kazakhstan,  Test  Site  (STS)  nuclear 
explosions  and  recorded  at  local  distances  by  the  Institute  for  Dynamics  of  the  Geosphere  (IDG). 
The  datasets  are  described  in  Stevens  et  al.  (2007)  and  has  excellent  examples  of  Rg  recorded  at 
local  distances.  The  test  sites  both  have  Rg  velocity  >  2  km/s  resulting  in  a  choice  for  the 
excitation  correction,  E(T)=-\.  The  MRg  estimates  for  all  18  events  were  converted  to  yield  using 
Eq.  9. 

The  comparison  of  the  true  and  estimated  yields  is  shown  in  Figure  6.  Because  of  the 
small  number  of  stations  recording  these  events  (typically  4-6  stations),  we  did  not  trim  the  mean 
for  these  results.  The  highest  percentage  error  between  the  estimated  and  true  yield  is  64%, 
which  suggests  a  smaller  F  factor  for  DTS  and  STS  that  what  was  determined  for  the  North 
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America  chemical  explosions  used  to  derive  Eq.  6.  For  12  of  the  18  events,  the  estimated  yield 
percentage  error  is  20%  or  less,  which  is  similar  to  Adushldn’s  (2001)  yield  estimates  for  Shagan 
River  Test  Site  nuclear  explosions  using  Rg  with  stations  corrections.  Considering  that  these 
DTS  and  STS  nuclear  explosions  have  yields  that  are  orders  of  magnitude  larger  and 
significantly  greater  depths  than  the  events  used  to  develop  Eqs.  6  and  9,  these  results  are  very 
promising. 

Over-buried  explosions  may  provide  difficulties  for  MRg  due  to  moment  decrease  and 
surface  wave  eigenfunction  depth  decay  for  Rg.  As  an  example,  Figure  5  shows  the  MRg  for  a  25 
ton  chemical  explosion  detonated  at  550  m  depth  (scaled  DOB  =  -1452  m/kt1'  )  as  part  of  the 
1997  STS  DOB  Experiment  (Myers  et  ah,  1999;  Patton  et  al.,  2005).  The  MRg  estimate  sits  well 
below  our  observed  data  population  of  events  with  sDOB  ranging  from  43  to  271  m/kt1/3  and 
results  in  a  yield  estimate  of  only  5  tons  for  this  25  ton  event.  However,  after  correcting  the 
observed  MRg  estimate  for  over-buried  moment  using  the  Denny  and  Johnson  (1991)  source 
model  and  surface  wave  depth  excitation  using  synthetic  seismograms,  the  corrected  MRg  is 
consistent  with  the  regression  equation  for  MRg  and  log  Y. 

Western  United  States.  We  have  experienced  difficulties  with  estimating  yields  for  Nevada 
Test  Site  (NTS)  explosions  using  the  MRg  method.  Analysis  of  near-source  recordings  (Brian 
Stump,  pers.  comm.  201 1)  of  NTS  explosions  such  as  Hunter’s  Trophy,  Mineral  Quarry,  and  the 
NPE  showed  little,  if  any,  observable  Rg  energy  at  nearby  stations,  perhaps  due  to  significant 
topographic  and  structural  scattering.  We  have  also  examined  recent  chemical  explosions 
conducted  at  NTS  (now  called  the  Nevada  National  Security  Site).  Because  1  s  Rg  was  difficult 
to  identify  at  distances  less  than  2  km,  we  modified  Eq.  5  for  application  on  3  Hz  Rg  recorded  at 
1-2  km  from  the  1000  kg  Source  Physics  Experiment  shot  2  (SPE-2;  Snelson  et  al.,  2012). 
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Results  of  the  MRg  methodology  varied  significantly  depending  on  path  effects,  as  we  estimated 
a  viable  yield  (892  kg)  for  three  stations  north  of  the  explosion,  while  yield  estimates  to  the 
northeast  ranged  from  486  kg  at  1.1  km  to  295  kg  at  1 .6  km.  These  results  demonstrate  that  there 
are  regions  where  the  scattering  and  attenuation  of  Rg  is  significant  and/or  highly  variable 
making  application  of  the  proposed  yield-MRg  methods  problematic  or  implausible.  With 
publication  of  this  paper  and  public  availability  of  MRg  estimation  software,  we  hope  other 
researchers  will  apply  the  formula  to  additional  explosion  datasets  and  inform  us  when  the 
method  fails  or  succeeds. 

5.  DATA  AND  RESOURCES 

The  software  developed  to  estimate  MRg  is  available  upon  request  from  the  authors.  Data 
used  in  this  project  is  derived  from  numerous  small-scale  explosion  experiments  conducted  by 
various  organizations.  Most  of  the  datasets  used  in  this  project,  which  are  older  than  2  years  from 
the  publication  date,  can  be  found  at  the  Incorporated  Research  Institutions  in  Seismology  (IRIS) 
Data  Management  Center  (DMC).  The  data  are  from  experiments  with  code  names  SPE,  FRE, 
NEDE,  HR,  BASE,  and  FDE.  For  data  that  are  not  located  in  the  IRIS  DMC  repository,  please 
contact  author  Bonner  for  details  on  how  to  obtain  those  data. 
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Figure  1.  Rg  attenuation  coefficients  as  a  function  of  period  for  diverse  lithologies  in  North 
America  as  well  as  data  for  Kazakhstan  from  Patton  et  al.  (2005).  Also  shown  are  the 
period-dependent  attenuation  terms  (Batt=0.0031  p=1.8)  for  the  Russell  (2006)  formula  as 
well  as  the  proposed  terms  (Batt=0.0031  p=1.96  and  Batt=0.0031  p=2.19)  for  T=1  s  Rg  based 
on  the  measured  attenuation  coefficients. 
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Figure  2.  Excitation  corrections  for  Ms  to  periods  as  low  as  0.5  s  for  five  different  crustal 
models  compared  to  the  excitation  correction  used  in  Russell  (2006). 
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Group  Velocity  (km/s) 


Figure  3.  Processing  example  showing  a)  filtered  waveforms,  b)  group  velocity  estimates 
for  Rg,  and  c)  magnitude  estimates  for  a  101.4  kg  (TNT  equivalent)  explosion  conducted  in 
Barre  granite  (Vermont,  USA). 
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Figure  4.  Final  MRg  estimates  (error  bars  show  ±la)  for  39  small  explosions  in  different 
lithologies.  The  results  are  plotted  against  the  TNT- equivalent  yields  for  the  explosions,  a) 
Results  using  three  excitation  corrections  (-1,  -1.5,  -2.5)  as  defined  by  Eq.  4.  b)  Results  using  a 
single  excitation  correction  (-1)  for  all  events. 
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Figure  5.  A  yield-MRg  relationship  based  on  regression  results  for  chemical  explosions.  The 
data  points  are  shown  as  solid  circles,  the  best  fit  regression  line  is  solid,  and  the  dashed  lines 
describe  the  95%  confidence  intervals.  The  open  circles  and  arrow  represent  results  from  an 
over-buried  Shagan  DOB  explosion  which  was  not  considered  in  the  regression.  The  arrow 
points  from  the  observed  Mr3  to  a  value  corrected  for  depth  of  burial  effects  due  to  moment 
and  Rg  excitation. 
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Figure  6.  Estimated  yields  from  MRg  (Y (%))  in  kt  versus  true  yields  for  8  Degelen  and  10 
Shagan  Test  Site  nuclear  explosions.  The  dashed  line  represents  a  1:1  relationship.  The 
vertical  error  bars  represent  an  F factor  of  2.25,  which  is  based  on  the  North  American 
chemical  explosion  datasets  used  to  develop  the yield-Mitg  relationship.  When  only  the  Shagan 
and  Degelen  nuclear  explosions  are  considered,  the  F factor  is  estimated  as  1.6. 
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